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Varying Cu-Ti hybridization near the Fermi energy in Cu2:TiSe2: Results from 

supercell calculations. 

T. Jeong and T. Jarlborg 
DPMC, Umversity of Geneva, 24 Quai Ernest- Ansermet, CH-1211 Geneva 4, Switzerland 

The properties of Cua;TiSe2 are studied by band structure calculation based on the density func- 
tional theory for supercells. The density-of-states (DOS) for x—0 has a sharply raising shoulder in 
the neighborhood of the Fermi energy, Ep, which can be favorable for spacial charge modulations. 
The Cu impurity adds electrons and brings the DOS shoulder below Ep- Hybridization makes the 
Ti-d DOS at Ep, the electron-phonon coupling and the Stoner factor very large. Strong pressure 
dependent properties are predicted from the calculations, since the DOS shoulder is pushed to higher 
energy at a reduced lattice constant. Effects of disorder are also expected to be important because 
of the rapidly varying DOS near Ep. 

PACS numbers: 71. 23.-15,71.45. Lr, 74.25. Jb,74.62.Fj 
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The layered TiSe2 compound is much studied because 
of the appearance of charge density waves (CDW) . The 
ground state is believed to be either a semimetal or a 
semiconductor with a small indirect gapfllSi. From 
photoemission experiments it is concluded that the CDW 
transition consists of a change in electronic structure 
from a small indirect gap into a state with a larger in- 
direct gap at a slightly different location in the Brillouin 
zonejl]. Recently Morosan et aZ.Q reported that a con- 
tinuous intercalation of TiSe2 with Cu is possible, where 
the Cu atoms enter between the TiSe2 layers to yield 
Cua;TiSe2, with solubility of Cu up to x « 0.11. The 
CDW transition is continuously suppressed, and a su- 
perconducting state appears near x = 0.04 with a max- 
imum Tc of 4.15 K at x = 0.08. A CDW supercon- 
ductivity phase diagram as a function of doping is devel- 
oped for Cua;TiSe2, analogous to the antiferromagnetism- 
superconductivity phase diagram found for the high tem- 
perature superconductors. 

In this work we investigate the electronic structure of 
Cua;TiSe2 for x=0 and x « 0.08 based on the density 
fimctional theory applied to supercell calculations. As 
will be shown, one of the main results is that Cu dop- 
ing introduces electronic charge carriers, leading to high- 
DOS properties, which partly confirms the interpretation 
of experimental results [5] . 

TiSe2 has the layered Cdl2 crystal structure. Planes 
of titanium atoms are surrounded by planes of selenium 
atoms, forming a sandwich structure, where the selenium 
atoms are arranged in octahedral coordination around 
the titanium atoms at z = ±i. The unit cells are ar- 
ranged so that the metal atoms are stacked directly above 
each other. The TiSe2 layers are well separated and 
bound by van der Waals interaction. As Cu atoms are 
added, they occupy positions within these open spaces 
between the TiSe2 layers. This results in a systematic 
expansion of the unit cell with increasing Cu content in 
Cua;TiSe2. The experimental lattice constants are a= 6.7 
a.u. and c=11.35 a.u. when x is near 0.08 [H]- 

The self-consistent Linear MufRn-Tin Orbital method 
[1, 01 is applied to a 48-site supercell in which one Cu 
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FIG. 1: The total density of states of TiSe2. The energy is 
relative to Ep. 



atom is intercalated between two Se layers. In these cal- 
culations we use a potential based on the local density 
approximation (LDA) 8], with no special on-site correc- 
tion due to correlation. The cell consists of CuTii2Se24 
and 11 empty spheres in the open space separating the 
TiSe2 layers. The repeat distance between Cu atoms is 
2 unit cells in the two plane directions and 3 cells along 
z. This is chosen to minimize the Cu-to-Cu interactions 
within a reasonably simple unit cell for a Cu concentra- 
tion X of 0.0833, which is close to the concentration for 
optimal Tc- The self-consistent iterations use initially 8 
fc-points in the irreducible Brillouin Zone, while they are 
terminated by iterations using 75 fc-points. The Wigner- 
Seitz sphere radii are 2.88 a.u. for Ti and 3.17 a.u. for Se 
respectively. The basis-set for the atoms includes d-states 
in the direct basis and f states as tail contribution. The 
direct basis for empty spheres includes p, with d states 
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FIG. 2: The total density of states of CuTii2Se24 at the lattice 
constant 6.7 a.u.. The energy is relative to Ef- 
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FIG. 4: The total density of states of CuTii2Se24 near Ef (at 
zero energy in each case) at three different lattice constants . 
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FIG. 3: The projected Ti-3d and Cu-3d parts of the total 
DOS shown in fig. 2. 



in the tails. The bands are calculated at the experimen- 
tal lattice constant, and at two smaller volumes, for the 
case with optimal Cu insertion. It should be recalled that 
LDA typically underestimates the lattice constant by 1-3 
percent for 3d metals. 

The supercell is repeated periodically and effects of 
random disorder, of distortion around the impurity, or of 
higher order periodicities among the Cu atoms are not 
included. It should be noted that any disorder (struc- 
tural or substitutional) within the supercell will intro- 



duce band broadening in general, since degenerate bands 
of the perfectly repeated supercell will be lifted. These 
effects of disorder are likely to be more important for 
this particular system, because of the very sharp DOS 
structures near Ep- 

The electron-phonon coupling A is calculated using the 
Rigid-Ion approximation and experimental information 
of the phonon moment (w^). As the Debye temperature 
for TiSc2 is not known to us, we use an average Deby e 
temperatures of NbSca (280 K f^) and TiSs (235 K ^). 
The modification of the pure phonon contribution to A 
due to the additions of Cu is probably not very impor- 
tant. Since the DOS is large on Ti, it is expected that the 
electronic contribution will be dominated by these sites. 
The pressure variation of the phonon moment is scaled 
through the calculated variation of the Bulk modulus. 
These procedures are not very precise, but they should 
be sufficient for getting insight of the contributions to A 

The Stoner factor, S, is calculated from the paramag- 
netic band results in the standard way, by assuming a 
finite rigid-band splitting due to exchange of all bands 
p^ . If there is only one band contributing to N{Ep), 
one can write S = N ■ I , where / is an exchange inte- 
gral for that band. The DOS is dominated by the Ti-d 
band and the Stoner factor is mainly because of Ti-d elec- 
trons. Hybridization with other bands, and the induced 
exchange splitting on other sites, is likely to reduce the 
effective Stoner enhancement, but spin-polarized calcu- 
lations would be needed to verify this. 

The DOS for the primitive cell of undoped TiSe2 is 
shown in Fig. 1. An independent calculation using the 
Full Potential Local Orbital (FPLO) [13] method, using 
a potential based on the gradient corrected density func- 
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tional |14| . gives a similar band structure, but with a 
slightly larger dispersion of a band above Ep, so that 
the first unoccupied DOS peak becomes less prominent. 
The band results indicate that the pure material is a 
semimetal where Ep is in a valley in the DOS, with a 
low value of N{Ef) of 10-20 states per Ryd per TiSe2. 
The uncertainty is due to the fact that Ep falls slightly 
on the right hand side of the minimum of the DOS (see 
Fig. 1), on the steep shoulder, where the DOS increases 
rapidly with energy. This is for the perfect structure 
without disorder or CDW. 

A simple rigid-band model for a potential (or charge) 
modulation in a material with an increasing DOS near 
Ep (as for pure TiSe2) can be shown to favor CDW. 
Shifts of the local DOS functions relative to Ep suggest 
that some regions of the material should get a high local 
N{Ep) because of the large DOS just above Ep. Oppo- 
sitely, other regions where the potential is repulsive will 
have Ep closer to the local "gap" with minimum DOS. A 
similar reasoning as for a Peierls transition show that the 
total kinetic energy will be lowered by such a splitting of 
local states. Namely, since the DOS is sharply raising at 
Ep and the number of states is conserved we find that 
an upward shift, 6u, of the DOS relative to Ep within 
the regions that loose electrons will be larger than the 
downward shift, dd, of the DOS within the region that 
gain electrons. If the average DOS to the left and right 
side are Nu and Nd, respectively, and An is the number of 
transfered electrons, we have that An = Su ■ = 6d ■ Nd- 
But the kinetic energy will be ^An{6d—Su), i.e., negative. 
Or equivalently, the total (kinetic) energy will be lower 
for the CDW state. If, by electron doping, Ep moves to 
higher energy on a linearly increasing DOS, then the rel- 
ative difference between Nu and Nd will be smaller, and 
the gain in kinetic energy and the tendency for CDW 
formation will diminish. 

Total energy contributions from the Coulomb interac- 
tions and hybridization are neglected here, but the model 
gives a hint that a sufficiently raising DOS near Ep is 
favorable for CDW. Nothing can be said about the peri- 
odicity of these waves from the model. However, from 
qualitative arguments it is expected that modulations 
with long wave lengths are more probable, since they 
have lower costs in terms of hybridization energy. It can 
also be noted that irregular regions of up- and down- 



shifted potentials, as for long-range disorder, will be fa- 
vored through the same type of energy gains as for CDW. 
Supercell calculations with some guidance for probable 
CDW modulations would be needed to solve this prob- 
lem. It is generally believed that Fermi surface nesting 
is not a determining factor for the k-vector of the CDW 

When Cu is inserted it acts as an electron dopant. Cu 
brings along d-states near Ep, which hybridize strongly 
with the Ti-d band, and the total DOS at Ep increases 
very much as can be seen in Fig. 2. From an inspection 
of the local Ti DOS functions it can be verified that the 
Ti close to the Cu site has the highest N{Ep) values. 
The total Cu-d and the sum of all Ti-d DOS functions 
are shown in Fig. 3. 

The DOS peak near Ep is high and narrow. This 
makes any calculation of N{Ep) dependent properties 
very sensitive to small shifts of Ep (non-stoichometry) 
and to band broadening (disorder effects). The energy 
derivative of N{Ep) is high, by shifting Ep 2 mRy to 
lower energies, N can be reduced by a factor of two. 

Nevertheless, by using the unshifted unbroadened band 
results we obtain a very large Stoner enhancement, S — 
1/(1 — ,5), see table I. Ferromagnetism is expected when 
S is larger than one. Apart from non-stoichometry and 
disorder effects it is possible that the impurity-like elec- 
tronic structure with a large DOS near the Cu site will 
overestimate the calculation of S. The reason is that 
the exchange splitting is assumed to be equal on all sites 
in the model, while in reality this may cost additional 
kinetic energy on sites with low DOS. A spin-polarized 
calculation would give a better answer, but it has not 
been attempted at this stage because of the very slow 
convergence for the supercell. 

The calculated values of A are reasonable and not as 
extreme as could be expected from the large DOS and 
S. The main contribution to A comes from the Ti atoms. 
The Ti-DOS is large, but there is only a moderate dipole 
scattering contribution to the electronic part of A, be- 
cause the total DOS is dominated by the d-character 
(Typically 98 percent of the Ti-DOS is of d-character). 

A surprising effect is the large variation of N{Ep) with 
pressure, P. Two calculations at smaller lattice con- 
stants show clearly that the narrow peak of the Cu-d 
state moves sensitively with volume. This state is hy- 
bridized with the Ti-d band so that also the Ti-d DOS 



TABLE L The calculated parameters of CuTii2Se24. The 
Debye temperature 6d is from experiment and scaled through 
the calculated volume dependence. The units of 7 is (mj/mol 
K^), and N(Ef) is in units of states per (Ryd-cell). 



TABLE 11: D-DOS at Ep on Cu and Ti near (Ticu) and far 
from (Ti) the the Cu impurity. The units of N(Ef) are states 
per (Ryd-atom) 



ao (a.u.) 


A 




S 


7 




N(Ef) 


6.7 


1.52 


20.2 


1.31 


11.13 


240 


1220 


6.67 


0.37 


0.11 


1.09 


4.75 


250 


960 


6.53 


0.21 


0.01 


0.49 


1.26 


290 


290 



ao (a.u.) 


Cu d N(Ef) 


Ticu-d N(Ef) 


Ti-d N(Ef) 


6.7 


22.1 


127. 


63.0 


6.67 


19. 


72. 


41. 


6.53 


0.85 


28.1 


1.0 



4 



at Ep will vary much with volume. A downward shift 
of this peak by about 5 niRy is found when the lattice 
constant is increased from 6.67 to 6.7 a.u.. If the lattice 
constant is reduced to 6.53 a.u. there will be an upward 
shift of the DOS edge to a location above Ep and the to- 
tal N{Ef) is not much larger than for pure TiSe2. The 
charges on Cu and Ti increases with P at the expense of 
charge on Se. The effect is that the Cu-d band becomes 
more filled for a uniform reduction of volume. The peak 
in the DOS 20 mRy below Ep for the smallest lattice 
constant, is mainly a Cu-d state. The Ti-d DOS is too 
far away in energy to hybridize with this state, and at 
Ep there will be less Cu-d to Ti-d hybridization. This 
will result in a reduction of N{Ep). An inspection of the 
local Ti-DOS for the smallest lattice constant reveals an 
important proximity effect near the Cu site. The Ti-d 
DOS at Ep is of the order 28 states/Ry/atom on each 
of the 8 Ti on the layers closest to Cu. But in the more 
distant Ti layer, where there are 4 Ti atoms, the local 
Ti-d DOS at Ep is only of the order 1-2 states/Ry per 
atom, which is like the case for pure TiSe2. The reason 
is that the local Ti-d DOS is pushed just above Ep for 
these sites. In the calculations for the two larger lattice 
constants the local Ti-d DOS values at Ep are compa- 
rable between all Ti sites, also within the most distant 
layer. This shows that the properties of Cu doped TiSe2 
will be very sensitive because of the exact position of the 
DOS shoulder relative to Ep. 

These calculations are made for perfect scaling of the 
supercell lattice. Thus no effects of possible rearrange- 
ment of the Cu occupations, different c/a ratio or disor- 
der are taken into account. If such effects are small, one 
can expect strong P-variations of the properties. Accord- 
ing to our calculations the system would pass from ferro- 
magnetism, towards exchange enhanced paramagnetism, 
through a region with superconductivity and large elec- 
tron phonon coupling, and finally to a fairly low-DOS 
material when the volume is reduced by about 8 percent 
only. However, this unusual sensitivity is caused by a 
sharp DOS shoulder which is passing through Ep, and 
therefore one can also expect large error bars. As men- 
tioned before, disorder and different Cu impurity distri- 
butions can be important and have a moderating effect 
on the P-dependence. 

More hybridization with the Ti-p and / bands is 
needed for a larger A. By reducing the lattice constant it 
is possible to increase the hybridization and the electronic 
contribution to A, but the effect is masked by the changes 
in total N{Ep), which is decreasing much at large pres- 
sure. In addition, fine structures in the peaked DOS add 
some noise to the evolution of N with pressure. 



An estimation of the superconducting Tc from the cal- 
culated electron-phonon coupling via the McMillan equa- 
tion gives reasonable values in line with observations, 
see Table I. However, some reduction is expected from 
spin-fluctuations, especially since the calculations indi- 
cate large exchange enhancement. The calculated elec- 
tronic specific heat constant, 7, includes the enhance- 
ment for electron-phonon coupling, but not from possi- 
ble spin fluctuations. At the largest lattice constant 7 is 
considerably larger than the measured value, about 4.5 
m J/mol K2 0. The calculated result for 7 at the in- 
termediate volume agrees best with experiment, but the 
electron-phonon coupling is rather low at the same vol- 
ume. Again, this is an indication that disorder in the 
real material will reduce the effective DOS at Ep and 
increase hybridization. 

In conclusion, we find that undoped TiSe2 is a low- 
DOS material with semi-metallic properties, but the DOS 
is rapidly increasing near Ep. The latter fact is probably 
important for the properies of pure and doped TiSe2 . For 
instance, the kinetic contribution to the total energy will 
be lowered by spacial modulations of the potential as pro- 
duced by a CDW or disorder. Electron doping, accom- 
panied by disorder, will put Ep closer to the broadened 
DOS peak. This should weaken the mechanism for CDW 
formation, while the material soon becomes a high-DOS 
material. The local Ti-d DOS at Ep for a supercell of Cu 
doped TiSe2 is increased very much through hybridiza- 
tion with the Cu-impurity and its d-band. The exact 
value of N{Ep) is uncertain because of the sharp struc- 
ture of the DOS near Ep, which is found both for pure 
TiSe2 and the supercell, and it depends probably also on 
Cu ordering. The estimation of A show that conventional 
superconductivity based on electron-phonon coupling is 
probable, with a Tc of the order of a few K. Despite the 
changes in electronic structure, from a low DOS semi- 
metal in pure TiSe2 to a high DOS material at high Cu 
doping, it is the Ti states that makes A large. Some re- 
sults, like the large Stoner factor, indicate that disorder 
in the real material is important to moderate the very 
large N{Ep). Effects of disorder will increase hybridiza- 
tion and can make A larger even if the total DOS is low- 
ered. The particular situation with Ep near a shoulder 
in the DOS, and electron doping via Cu- impurities, turns 
out to give strong volume dependence of N{Ep) in ad- 
dition to the dependence on x. Further work is needed 
to sort out the precise effects of disorder, but it can al- 
ready be predicted that not only variations of x, but also 
P-variations, will be determining for the properties of 
Cua;TiSe2. 
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